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ABSTRACT 

We present a candidate orbital period for the low mass X-ray binary XB 1832-330 in the globular 
cluster NGC 6652 using a 6.5 hour Gemini South observation of the optical counterpart of the system. 
Light curves in g' and r' for two LMXBs in the cluster, sources A and B in previous literature, were 
extracted and analyzed for periodicity using the ISIS image subtraction package. A clear sinusoidal 
modulation is evident in both of A's curves, of amplitude ~0.11 magnitudes in g' and ~0.065 mag- 
nitudes in r', while B's curves exhibit rapid flickering, of amplitude ^1 magnitude in g' and ~0.5 
magnitudes in r'. A Lomb-Scargle test revealed a 2.15 hour periodic variation in the magnitude of A 
with a false alarm probability less than 10 -11 , and no significant periodicity in the light curve for B. 
Though it is possible saturated stars in the vicinity of our sources partially contaminated our signal, 
the identification of A's binary period is nonetheless robust. 

Subject headings: binaries : X-rays — stars: neutron — accretion — globular clusters (individual): 
NGC 6652 



1. INTRODUCTION 

Globular clusters are efficient factories for the produc- 
tion of low-mass X-ray binaries (LMXBs), due to their 
high cen tral densities and increased rates of s tellar inter- 
actions (IClarkl [19751 : IVerbunt fc Llwir][2006Tl . The dy- 
namical formation mechanisms in globular clusters pro- 
duce LMXBs with different characteristics than those in 
the rest of the galaxy, particularly in their overproduc- 
tion of ultracompact (P <80 minutes) LMXB s requiring a 
degenerate donor star (jDeutsch et al.l [20001 ). Numerous 
LMXBs have been detected in globular clusters around 
other galaxies (e.g. iSarazin et alJ 120011 : lAngelini et alJ 
l2001f ). allowing ident ification of a metalli city dependence 
in their formation (|Kundu et alJ I2002D . better statis- 
tics o f their dependence on cluster structural param- 
eters ([Jordan et al.l [2007T). and high-quality X-ray lu - 
minosity functions (|Kim et alJ 120091 : iZhang et all 120111 ) . 
However, local globular cluster LMXBs are generally the 
only objects we can study in detail, obtaining key pa- 
rameters such as binary orbital peri ods for use in pop- 
ulation synthesis c alculations (e.g. [Fragos et alJ 120081 : 
llvanova et all 120081 ). As of 2011, 15 bright (L x , pe ak > 
10 36 ergs/s) LMXBs in galactic globular clusters are 
known; of eleven known periods, five are ultracompact, 
while two require (P or b_ >15 hours) evolved donor stars 
(IVerbunt fc Lewinll2006l: iDieball et apl2005l:lZurek et al ' 



20091: lAltamirano et ail 120081 . I20iq Istrohmaver et al 



20101 ). Identification of these orbital periods has gener- 



ally required Hubble Space Telescope (HST) imaging, or 
X-ray studies, w ith the exception of t he optically bright 
AC 211 in M15 ([Ilovaiskv et al.|[i"993l ). 

X-ray emission from the globular clu ster NGC 6652 
was detected by HEAO-1 in 1977-78 (|Hertz fc Woodl 
[l985h . at i x (2-10 keV) - 10 35 ergs/s, and then at L x ~ 
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10 36 ergs/s by ROSAT in 1990 ([Predehl et al.lll99l . 
This LMXB has sinc e been o bserved by BeppoSAX 
din 't Zand et~aT1 fl998l) ASC A (|Mukai fc Smald l200l . 
and XMM ([Sidoli et al.ll2008f ) at L x ~ 10 36 ergs/s. Since 
1999 it has been monitored by RX TE's PCA instrument 
during the bulge scan program (|Swank fc Markwardtl 
l200ir FI. finding it roughly constant (Lx ~ 10 36 ergs/s) 
u ntil 2011, during whi ch it declined by a factor of ~3. 

iDeutsch et al.l (|1998l) identified a candidate UV-bright, 
variable optical counterpart (their star 49) to the 
NGC 6652 LMXB a t 2.3ct from the ROSAT position. 
IDeutsch et alJ (|2000l ) found a candidate 43.6 minute pe- 
riod in three orbits of HST V and / imaging, with a 
Fourier peak at 99.5% confidence. However, only two 
HST orbits (each including ^45 minutes of on-source 
time) are consistent with this period, while the third 
orbit shows s trong flickering (see their Fig. 1), and 
iHeinke et al.l (120011 ) argu ed the period was not convinc- 
ing. IDeutsch et alJ (|2000l ) note that the optical faintness 
of star 49 indicates a small accretion disk, and thus also 
suggests an ultracompact nature. 

Chandra high-resolution X-ray imaging of NGC 6652 
has revealed several X-ray sources in the cluster, here- 
after A through G in orde r of brightness (jHeinke et alJ 
12001 iCoomber et al.ll201ll ). The bright LMXB A (XB 
1832-330) was identified by IHeinke et all ([20011 ) with 
a different blue, variable star, showing tentative evi- 
dence for a sinusoidal period of either 0.92 or 2.22 hours. 
Deutsch's star 49 is the optical counterpart of B, a lower- 
luminosity (Lx ~ 4 x 10 33 ergs/s) X-ray source. A 5 
kilosecond Chandra exposure revealed rapid flaring from 
B on timescales down to 100 seconds, from Lx < 2 x 10 33 
ergs/s up to Lx ~ 10 35 ergs/s. B's long-term Lx has ap- 
peared relatively constant, as it has been detected with 
ROSAT in 1 994, and Chandra in 2000, 2008, and 2011 at 
similar L x s (jCoomber et al.ll2011i IStacev et al.ll2011[ ). 

A and B are somewhat unusual among globular clus- 
ter LMXBs in being located well outside the core of their 

4 http://lheawww.gsfc.nasa.gov/users/craigm/galscan/main.html 
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globular cluster (|Verbunt fc Lewinll2006f) , in less crowded 
regions potentially observable from the ground. This mo- 
tivated us to search for optical periodicities in both A and 
B using Gemini-South. 

2. DATA REDUCTION 

We observed NGC 6652 on 2011 May 2 for 6.5 hours 
with Gemini's GMOS-S CCD detector. Each exposure 
was 75 s in duration, and 172 images were taken alter- 
nately in the g' and r' filters. Raw CCD data were pre- 
pared by the Gemini IRAF0 package task GPREPARE, 
and flat-field and bias corrections, as well as gain multi- 
plication, were performed by the task GIREDUCE using 
Gemini-supplied calibration information^. The GMOS-S 
instrument contains three CCD chips, but only the im- 
age data from the central chip, containing the majority 
of the cluster, was used. 

Photometry was performed on the images using the 
image s ubtraction package ISISlL developed by Alard & 
Lupton (jAlard fc Luptonlll99alAlardll2000h to deal with 
very crowded fields exhibiting spatially varying seeing 
and background levels. ISIS works by transforming im- 
ages to a common seeing, subtracting these from one an- 
other, and performing photometry on the more sparsely 
populated subtracted images. To minimize PSF residuals 
in the subtracted image, ISIS uses a linear least squares 
fitting method to optimize a solution for a spatially- 
varying kernel which can be convolved with a reference 
image to recreate the seeing of individual images. The 
basis functions of this kernel are taken to be the products 
of polynomials and Gaussian distributions. 

The steps ISIS performs are as follows: first, all images 
are transformed to a common template coordinate grid 
using a polynomial astrometric transformation. Next, 
several of the images with the best seeing are transformed 
to the same seeing and stacked to form a reference im- 
age. A best-fit kernel is then computed for each image, 
which is convolved with the reference to match the ref- 
erence image to the seeing of each image. Each image is 
then subtracted from the convolved reference image. The 
subtracted images are stacked to form a median image; 
bright spots in this image indicate either variable stars 
or saturated stars, for which the PSF matching cannot 
be done properly, and which therefore generate residuals. 
PSF-fitting photometry is performed on each subtracted 
image to generate a light curve for variable objects. 

Using a trimmed version of the image (151 x 301 pixels) 
with ISIS avoided undesirable effects of the many satu- 
rated foreground giant stars in the full (1024 x 2304 pix- 
els) image. Since these trimmed images were a manage- 
able size, they were processed in a single section (sub_x 
= sub_y = l). 

Because the coordinates of the two stars of interest, 
A and B, were already known, ISIS parameters which 
optimized the variation signal for these sources in the 
median stacked subtracted image (abs . fits) and also 
minimized contamination by PSF residuals of nearby 
saturated stars (see section 4) were chosen. Each im- 
age was transformed to a common coordinate grid using 
two-dimensional polynomial astrometric re-mapping of 

5 http://iraf.noao.edu/ 
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degree 1 (DEGREE = 1). 

In the g' filter, a reference image for subtraction was 
created using the four images with the best seeing, and 
in the r' filter, eleven best-seeing images were used. Im- 
ages were processed using 10 stamps in each direction, 
each of radius 15 pixels (nstamps_x = nstamps_y = 
10 , half_stamp_size=15), a radius for the convolu- 
tion kernel of 9 pixels (half_mesh_size = 9), third 
degree variations in the background level (deg_bg=3) 
and a second degree spatial variation of the kernel 
(deg_spatial=2).The saturation level was left at 1 000 
000, as lowering it seemed to contaminate our lightcurves 
more severely. After image subtraction, median stack- 
ing, and variable detection, the ISIS-derived coordinates 
(in phot, data) for the centres of A and B were ad- 
justed manually so as to be more aligned with the sources 
prior to photometry. Photometry was then performed 
on the subtracted images using a fitting radius of 6 pix- 
els (radphot = 6.0). An exploration of the results of 
varying other photometry parameters led us to conclude 
that the default values supplied by ISIS were sufficient 
for our purposes. 

3. LIGHT CURVE EXTRACTION & CALIBRATION 

The lightcurves from the ISIS photometry for A and 
B appear in Figure [TJ Calibration was performed 
in two steps, a c cordi ng to the method employed by 
iMocheiska et al.l (|2001h : first, the magnitudes of A and B 
in a template image were compute d using the DAOphot 
package ALLSTAR (|Stetsonlll987| ). An aperture correc- 
tion was applied to these magnitudes by comparing the 
flux admitted through apertures equal to the PSF and 2 
pixels larger than the PSF for stars in the vicinity of A 
and B (central region of the image in Figure 0]). 

The template image in both filters was chosen to be 
one with superior seeing. The template magnitudes, 
m tp i , for each source were converted into template im- 
age counts, c tp i, using the ALLSTAR zero point of 25.0 
magnitudes and the exposure time of 75.0 s using the 
relation: 



m tpl = 25.0- 2.5 log(^) 



(1) 



Once the counts for A and B in the template image were 
obtained, it was possible to determine their counts in 
the stacked reference image used for subtraction, c rG f , by 
adding the counts on the template image and the delta 
count value ISIS reported for the subtracted template 
image, Ac tp i: c ro f = ct p i + Ac tp i. To convert the light 
curve point by point into magnitudes, the count values 
for A and B were computed for each image by subtracting 
the image's ISIS delta count values from the reference 
image count values: 



Cref - Ac, 



and inserting the results into equation Q] Errors in mag- 
nitude were computed from the ISIS flux errors using the 
relation 

6m = 2.5 log(l ± 6 JBH^L) (3 ) 
counts 

where the upper error bound corresponds to the addition, 
and the lower, to the subtraction. 

The second part of the calibration involved attempting 
to compute the 'true' apparent magnitudes using com- 
parison with photometric standards. Two standard field 
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Fig. 1. — Top left: ISIS g' light curve for A, with two fits which 
are combinations of a sinusoid, a linear function, and a constant 
offset. Fitting parameters were computed by QDP for the blue 
line shown. The (similar) red line is identical save for its period, 
which is the output from a Lomb-Scargle test. A secular trend 
of decreasing brightness is also apparent. Bottom left: The same 
as above, but for A in the r' filter. The variation is less obvious, 
but as discussed in the text, the two periods represented by the 
red and blue lines agree quite closely to those apparent in the g' 
filter. Top right: B's g' ISIS light curve. No periodic variation is 
evident. Bottom right: B's r' ISIS light curve. Again, no periodic 
variation is apparent. Errors shown here represent relative errors 
computed by ISIS only; they do not reflect errors associated with 
the conversion into magnitudes. 

images, of E8-a and 160100-600000, taken with the same 
CCD chip on the morning of May 2, 2011, were down- 
loaded from the Gemini ar chive. The Souther n u'g'r'i'z' 
standard star catalog from lSmith et al.l (|2007f ) was used. 
A total of 16 standard stars in g' and 20 in r' from across 
both fields were compared to derive calibration factors 
of 3.12±0.08 mags in g' and 3.25±0.04 mags in r' after 
weighted averages were taken. It should be noted that 
the factors derived for each field separately varied on the 
order of 0.1 magnitude for both filters. 

4. LIGHT CURVES & CONTAMINATION EFFECTS 

While Figure _] shows no periodic fluctuation in B's 
light curve, it indicates a clear periodic variation for A 
in both filters. When A's light curve is compared to the 
ISIS light curves for other stars, however, it becomes ap- 
parent that a possible contamination exists. In Figure _3 
the g' light curves for several stars, including A, and the 
inversion of the time variation of the PSF as computed 
by the IRAF task PSFMEASURE are shown. The PSF 
has been inverted here for visual aid in identifying curve 
similarities. Most of the saturated stars in the image 
displayed light curves that closely resembled the inverse 
PSF, much like the curve of Star 57 in the figure. This is 
understandable, as the PSF fitting routine in ISIS can- 
not properly handle saturated stars, and thus significant 

8 http://www-star.fnal.gov/Southern_ugriz/www/Fieldindex.html 



ringlike residuals - due to incorrect PSF weighting - ap- 
pear in the subtracted images. When the PSF is large, 
the computed flux for saturated stars is too small. 

Light curves of several unsaturated stars in the vicin- 
ity of saturated stars also appear contaminated by the 
PSF resi duals of their neighbou rs; see Stars 37 and 91 in 
Figure _3 lHartman et alj (|2004f ) noted a similar spurious 
variability introduced in the light curves of stars near 
saturated sources. Not all unsaturated stars near satu- 
rated stars in our images exhibit this effect, however; the 
curve for Star 63 (see Figure [2]), which is adjacent to sat- 
urated star 57 (see Figure _|, displays little similarity to 
the PSF curve. Looking at the curve for A, it seems plau- 
sible that while the PSF has impacted the result to some 
degree, a true periodic signal exists. Marked dissimilari- 
ties between the locations of the peaks in the PSF curve 
and source A's curve corroborate this assertion. Results 
in the r' filter are similar. 

It is evident from Figure _] that the comparison stars 
in Figure _0 are all brighter than A. It is desirable to 
compare A to a faint star in the vicinity of a saturated 
star; however, ISIS photometry can be performed only 
on variable stars, and few of the faint stars in the image 
were detected as variable by ISIS. Two faint stars that 
were picked up by ISIS on the subtracted images, Stars 
87 and 88, appear in Figure _3 While these are farther 
from saturated stars than A is, their smaller amplitude 
variations (compared to A) suggests that A's signal is 
not spurious. Indeed, Star 91's lightcurve (also shown in 
Figure [3]) is also of smaller amplitude than A's, despite 
its proximity to the saturated Star 84. Finally, all of 
the comparison lightcurves in Figures _0 and [3] possess 
greater scatter than A's relatively smooth curve. We 
therefore concluded that the sinusoidal variation in A 
is a reflection of the binary orbit - possibly distorted 
somewhat by the saturated stars' PSF residuals - and 
proceeded to compute a candidate period for the motion. 
Figure _] presents the star field of interest, along with an 
enlarged image displaying stars used for comparison in 
Figure O 

5. ANALYSIS OF PERIOD 

To determine the period of A's variation, a Lomb- 
Scargle Periodogram (iLomrj[T976t IScargldfl982T ) was cal- 
culated for data in each filter. For the g' data, a clear 
peak appeared at a period of 2.148 ± 0.002 hours with 
a false alarm probability (FAP), or probability that the 
data was drawn randomly from a Gaussian distribution, 
of 8.2 x 10~ 12 . To determine the error in the period, 
a Monte Carlo method was employed. First, the errors 
in each data point were multiplied by a factor randomly 
drawn from a normal distribution about with a stan- 
dard deviation of 1.0. The result of this multiplication 
was then added to each data point, and the Lomb Scar- 
gle test was applied to the new data set. This process 
was repeated 10 4 times using an oversampling param- 
eter of 2048, sufficiently high to adequately sample the 
frequency space. The results for the peak period for each 
of the 10 4 trials were plotted and fit to a pseudo-Gaussian 
distribution; the values on both sides of the central peak 
which were equivalent to the la mark in terms of outlying 
percentage of data (~15.8% lying outside these limits on 
either side) were taken to represent the upper and lower 
error limits. Figure [5] displays the results of our Monte 
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Fig. 2. — Lightcurves for several stars in the g' filter, rescaled 
and shifted. To emphasize curve resemblances, the inverse of the 
PSF curve is plotted at top. Star 57, a typical example of a sat- 
urated star, has a curve which closely mirrors the inverted PSF. 
Unsaturated stars in the vicinity of saturated stars can also be af- 
fected, as seen in Star 37's light curve and Star 91 's light curve, 
which follows the PSF closely. Some unsaturated stars near bright 
stars do not show this effect, such as Star 63 and, to some extent, 
Source A. 

Carlo trials. 

The above analysis was repeated for the r' filter data 
to yield a period of 2.149±0.004 hours, with an FAP of 
2.7 x 1CT 9 . 

The Lomb-Scargle periodograms for several stars in the 
field are shown in Figure [6] The two highest peaks cor- 
respond to A in the g' and r' filters, respectively. Star 
91, itself unsaturated, is located above the saturated star 
adjacent to A (see Figure [4}, and thus likely suffers from 
the same PSF contamination that A experienced, leading 
to the peak in its periodogram somewhere between that 
of A and the PSF. Lack of complete coincidence with 
the PSF period peak could be attributed to the photom- 
etry of star 91 picking up the periodic signal from A, 
though the distance between the stars presents an ob- 
stacle to this interpretation. In any case, the lightcurves 
in Figure [2] present convincing evidence that any peri- 
odicity detected in star 91 is spurious, and an effect of 
contamination. Peaks for other stars we examined were 
less significant than those shown. For visual aid, sev- 
eral FAP levels are plotted alongside the periodograms 
in Figure [BJ To determine the Lomb power these FAPs 
correspond to for plotting purposes, the relation 



Pi 



FAP 



z\M 



(4) 



was used, where Pfap is the false alarm probability, z is 
the Lomb power, and M is e quivalent to the n umber of 
data points for our purposes (|Press et alJll992T ). 

No significant peaks appeared when a Lomb-Scargle 
test was performed on B's light curves, as shown by the 
right panel in Figure [6] 



Fig. 3. — Lightcurves for several stars in the g' filter shifted, 
but not rescaled. Stars 87, 88, and A appear similar in brightness 
in Figure [4] yet the amplitude of the variations in A's lightcurve 
are much more significant. Star 91 is much brighter than A, as 
seen in Figure [4] and yet exhibits smaller amplitude variations, as 
well. The scatter in the curves for stars 87, 88, and 91 is also much 
more significant than that for A; this and the amplitude contrast 
evident above lend confidence to the interpretation of A's variation 
as a real signal. 
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Fig. 4. — Left: Portion of the image from the central CCD chip 
of Gemini's GMOS-S camera. The rectangle marks the region at 
right. Right: Closer view of the region of interest. Stars whose 
light curves appear in Figures [2] and [6] are labeled, as well as A and 
B. The bright star between A and Star 91 is saturated, as are the 
stars to the left of Star 37. North is toward the top of the image. 
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Fig. 5. — Left: Distribution of periods of 10 4 resampling trials for 
Lomb Scargle period of A using g' filter data. Errors for each data 
point were multiplied by a factor randomly drawn from a normal 
distribution about with a standard deviation of 1.0, then added 
to each data point. A Lomb Scargle test, using an oversampling 
parameter of 2048, was performed on each resulting data set. The 
horizontal line shows an approximate 1 a error designation, such 
that ~15.8% of the trials lie on either side of the bounds. Right: 
Same as detailed above, but for the r' filter data. 

A separate method ol fitting the light curves was pur- 
sued to obtain an independent estimate of the periodic- 
ity. The plotting program QDF0was used to fit A's light 
curves to the functional form 



27T 

Asin(— + C)) +Dx + E. 
B 



(5) 



Due to an initially large reduced x 2 value obtained upon 
fitting, the errors were increased until the % 2 was com- 
parable to the degrees of freedom. Fitting the data set 
with seven times the initial errors yielded parameters A 
= 0.119±0.008; B = 2.15±0.02; C = 7.01±0.02; D = 
0.031±0.003; and E = 19.40±0.02. This fit is plotted in 
Figure [1] in blue. The same fit with the period changed 
to the result of the Lomb-Scargle test for comparison is 
shown in red. 

An identical fitting procedure was applied to the r' 
filter data, with the results A = 0.065±0.009; B = 
2.16;°;^; C=0.5±0.2; D = 0.016±0.004; and E = 
19.69±0.03. Errors here are reported at the la level. 
Again, the fit (blue), along with an altered version con- 
taining the Lomb-Scargle period (red), appears in Figure 
ffl 

We are confident in our estimate of the candidate pe- 
riod due to the agreement of the Lomb-Scargle period es- 
timations in both filters with one another and the agree- 
ment of the QDP fit results with the Lomb-Scargle esti- 
mates within la. 



Fig. 6. — Left: Lomb-Scargle periodogram displaying the signif- 
icance of the periodicities in the light curves of several stars. The 
highest peaks belong to A in the g' and r' niters, respectively. Also 
shown is Star 91, an unsaturated star in the vicinity of A and satu- 
rated star 84; periodicity in its light curve suggests a contaminating 
effect of PSF residuals near saturated stars, corroborated by the 
shape of its lightcurve in Figure[2] A periodicity with FAP ~ 10~ 4 
appears in the PSF variation, and the curves of saturated stars, of 
which star 57 is a typical example, display the same peak. The four 
dotted horizontal lines mark confidence levels; the FAP for each is 
indicated. Right: Lomb-Scargle periodograms for B in both filters, 
accompanied by periodograms for the PSF and A for reference. No 
significant periodicity is seen in B's optical light curves. 

6. COLOR MAGNITUDE DIAGRAM 

The DAOphot package ALLSTAR was used to per- 
form PSF fitting photometry on the best-seeing images. 
Color-magnitude diagrams were constructed using im- 
ages taken consecutively in g' and r', separated by ~four 
minutes. Due to severe crowding by saturated stars (see 
Figure [4j , we could obtain usable photometry for A and 
B in only a few images. Figure [7] shows our colour- 
magnitude diagram with the tightest main sequence that 
shows A and B. A caveat is that the reduced % 2 values for 
A and B, indicative of goodness of PSF fit in DAOphot, 
were markedly higher (7.0, 5.1, 7.1 and 3.8 for A and B 
in r' and g', respectively) than for most main sequence 
stars (typically 1.0), indicating that the photometry for 
A and B is less reliable than that for most main sequence 
stars. 

The position of the main sequen ce turn-off is compa- 
rable to that in lHeinke et alj (|2001[ ). B again lies on the 
main sequence, while A is bluer than it, similar to results 
in those authors' HST color-magnitude diagram. 

7. DISCUSSION 



The possible contamination of A's source lightcurve by 
saturated stars forces us to exercise caution in our attri- 
bution of a period. However, several arguments point to 
http://heasarc.gsfc.nasa.gov/docs/software/ftools/others/qdp/qdp.li|jnl 2.15 hour period being real. The Lomb power peak 
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Fig. 7. — Color-Magnitude Diagram constructed using consecu- 
tive images in r' and g' with good seeing. All stars within 66.5" of 
the cluster centre were included in the image, and the breadth 
of the main sequence indicates the typical errors. Due to large \ 2 
residuals given by ALLSTAR for A and B, their magnitudes should 
be regarded as suffering larger errors. 

for the PSF modulations and other saturated stars was 
~1.9 hours, significantly different from A's Lomb peak. 
Figure [2] shows that the shape of the lightcurve due to 
the PSF modulation is significantly different from that 
of A, and Figure [3] displays A's larger amplitude and 
smoother variation with respect to stars similarly situ- 
ated. A sinusoidal shape is a very good description of 
A's lightcurve, as seen in Figure [1] Finally, the 2.15 
hour period is seen clearly in both g' and r' , and is con- 
si stent with the l onger of the two candidate periods for A 
in lHeinke etaLl (|2001l) . iHeinke et all (|200lD 's paper rep- 
resents an independent dataset confirming our own ob- 
servations; furthermore, since their data was taken with 
HST, for which neither seeing nor saturated stars were an 
issue, the agreement lends weight to our interpretation 
of the 2.15 hour period as a real signal. A repeat ob- 
servation, using a shorter frame time to prevent nearby 
stars from saturating, could give final confirmation of 
this period. 

A 2.15 hour period and A's persistent X-ray luminosity 
of ~ 10 36 ergs /s for 20 years suggest a system similar 
to GS 1821-26 (jHomer et alJll99ctlMeshchervakoy et al.l 
[20lB and a donor star of order 0.2 M g (cf. lDelovdl2008D . 

A sinusoidal modulation of ~ 0.1 magnitude could be 
due to heating of one face of the donor, or ellipsoidal 
modulations due to the varying donor aspect, as o f- 
ten seen for cataclysmic variables (|Edmonds et al.lfeOOSj ) , 
which would imply a 4.3 hour period. We rule out el- 
lipsoidal variations on two counts, however: firstly, the 
larger amplitude of the g' vs. r' modulations of A im- 
ply that the surface with changing visibility is quite hot, 
like the face of the donor heated by the disk. Further, 



a 0.2 M Q star should be located near the bottom of the 
main sequence. A's position (see Figure [7]) suggests that 
there must be substantial contributions to the luminosity 
from the bluer accretion disk and the heated face of the 
secondary. The relative variation of flux, and therefore 
the magnitude variation, i s much smaller than it would 
be were no disk visible. iWang fc Chakrabartvi (|2004f ) 
outline a case with a similar variation amplitude corre- 
sponding to sinusoidal variations. The small amplitude 
of A's variations may be explained by emission of the 
disk overwhelming the sinusoidally modulated visibility 
of the hot side of the donor star. We therefore think the 
identification of 2.15 hours as the binary orbital period 
is solid. 

A's position in the g' , g'—r' color magnitude diagram is 
slightly lower than the main sequence turnoff , in contrast 
with its posi tion slightly above t he turnoff in a V, V — I 
diagram by IHeinke et alJ l|2001f l. We expect that the 
majority of the light in this system is produced by the 
hot accretion disk, so the X-ray decay observed in 2011 
(see Introduction) may explain the decrease in the optical 
brightness. 

B's lightcurve (Figure [T]) appears completely chaotic. 
While this suggests that ISIS may not be producing reli- 
able results, this flickering is more likely to be intrin- 
sic to the system. The flickering is larger in g' (~1 
mag) than r' (~0.5 mag), which is expected if it relates 
to X-ray reprocessing, but not for saturation effects (as 
nearby giants are more saturated in r'). The flickering 
is not unlike the variations seen bv lDeutsch et al.l (|2000h 
in B's HST lightcurve, both in amplitude (~1 mag in 
V and I) and timescales (5-20 minutes). Finally, B's 
Chandra X-ray lightcurve shows order-of- magnitude flar- 
ing with typical timesca les of 5 minutes (|Coomber et al.1 
120111 : IStacev et aLll2011| ). Considering these arguments, 
we think that the ISIS photometry for B may be ac- 
curately reflecting the optical variation of the system, 
making it an intriguing target for future simultaneous 
X-ray /optical studies. If B's photometry presented here 
is reliable, it strongly indicates that the 43 minute candi- 
date orbital period of iDeutsch et al.l (|2000f ) is spurious. 

The optical color of B remains a significant puzzle. 
Since B's flickering occurs on the same timescales as 
our exposure lengths, it's difficult to accurately measure 
its color; however, the g' and r' measurements suggest 
a color as red as, or redder than, the main sequence. 
The flaring requires some component of the system-the 
disk, the companion, or both-to be strongly heated by 
X-rays, which should make it blue. Thus, B's color sug- 
gests a donor redder than the main sequence, perhaps a 
subgiant or "re d straggler" star with an unusual evolu- 
tionary history flAlbrow et al J 120011 iFerraro et all 120011 : 
iMathieu et alJl2003| ). 

8. CONCLUSION 

We have identified a clear 2.15 hour sinusoidal mod- 
ulation in the g' and r' lightcurves from the LMXB 
A in NGC 6652. Although contamination of the ISIS 
lightcurves due to signals from saturated stars is possi- 
ble, we are confident the sinusoidal signal is robust and 
likely represents the binary orbital period for reasons out- 
lined above. Further Gemini imaging, with shorter ex- 
posure times to avoid saturating nearby stars, could fi- 
nalize our result. We note that this period, if confirmed, 
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will be only the second X-ray binary in a globular clus- 
ter to have its orbital period measured using ground- 
based tel_e^car2es J _after_Jthe much brighter) AC 211 in 
M15 (|Ilovaiskv et alJll993h . 

The low-luminosity LMXB, B, in NGC 6652 shows a 
highly variable lightcurve with strong ~1 magnitude vari- 
ations on timescales of 5-20 minutes. This optical flicker- 
ing (if real) is probably driven by the strong, rapid X-ray 
flaring seen from B. 
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